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I. INTRODUCTION 

Directed polymers in a quenched random potential have been the subject of intense investigations during the past 
two decades (see e.g. PHI])- In the one-dimensional case we deal with an elastic string directed along the r-axis within 
an interval [0, t]. Randomness enters the problem through a disorder potential V[0(r), r], which competes against the 
clastic energy. The problem is defined by the Hamiltonian 



I' 1 

H\6(t),V] = 



f dr{\[d T <l>(r)] 2 + V[^r),r}}; (1) 



where the disorder potential V[</>, r] is Gaussian distributed with a zero mean V(4>,t) — and the ^-correlations: 



V(<t>, t)V{<P, t 1 ) = u5{t - T')8{(j> - <j>') (2) 

Here the parameter u describes the strength of the disorder. Note that such system is equivalent to the problem of 
the Kardar-Parisi-Zang (KPZ) equation describing the growth in time of an interface in the presence of noise [7j. 

In what follows we consider the problem in which the polymer is fixed at the origin, 4>(0) = and it is free at r = t. 
In other words, for a given realization of the random potential V the partition function of the considered system is: 

Z= I dxZ(x) = exp{-/?F} (3) 



where 

Z(x) = / Xty(r) e~ m ^ ] (4) 

J0(O)=O 

is the partition function of the system with the fixed boundary conditions, (f)(0) = and <j>(t) = x and where F is 
the total free energy. Besides the usual extensive part fot (where /o is the linear free energy density), the total free 
energy F of such system is known to contain the disorder dependent fluctuating contribution F. In the limit of large 
t the typical value of the free energy fluctuations scales with t as F oc t 1 ^ 3 (see e.g. 3-6]). In other words, the total 
free energy of the system can be represented as 

F = f t + ct 1 / 3 f (5) 

where c is a non-universal parameter, which depends on the temperature and the strength of disorder, and / is 
the random quantity which in the thermodynamic limit t — > oo is described by a non-trivial universal distribution 
function P(f). Note that according to eqs.((3|)-((5]), the trivial self-averaging contribution fot to the free energy can be 
eliminated by a simple redefinition of the partition function: 

Z = ex p{-pf t}Z (6) 
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so that 

Z = cxp{-A/} (7) 

where 

A = /3ci 1/3 (8) 

For the similar problem with the zero boundary conditions, 4>(0) = <f>(t) = 0, the corresponding distribution function 
was proved to be described by the Gaussian Unitary Ensemble (GUE) Tracy- Widom distribution |8l4ll|. In the course 
of this proof rather efficient Bethe ansatz replica technique has been developed [13, EJ. In particular, in terms of 
this technique the corresponding multi-point free energy distribution functions have been derived [l2j]. Recently, the 
free energy distribution function for the directed polymer problem with the free boundary conditions, eqs.([T]) -((][]), has 
been obtained [13j . It was shown that the function P(f) is the Gaussian Orthogonal Ensemble (GOE) Tracy- Widom 
distribution. In this paper I would like to present sufficiently simple alternative way of derivation of the same result 
which does not require rather complicated technique of the Fredholm Pfaffian described in [l3[ . 

Let us introduce the function 

y* oo 

W(f) = / dfP(f) (9) 
Jf 

which gives the probability that the random free energy is bigger that a given value /. It will be shown that in the 
thermodynamic limit, t — > oo, this function is equal to the Fredholm determinant 

W(f) = det(l - £_,) = Ft(-f) (10) 

with the kernel 

K-f{u,(J) = Ai(w + «'-/) ; (w,w' > 0) (11) 
which is the GOE Tracy- Widom distribution 0, Ell- Explicitly, 



Fx(s) = exp 



i ^ + 00 -i r + oc 
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(12) 



where q(£) is the solution of the Panleve II differential equation, q"(£) = + 2q 3 (£), with the boundary condition 

9 (£-> +0O ) = Ai(0- 

It should be noted that the present paper is rather technical. The main message of this work is not the final result 
itself (which is well known anyway) but the presentation of the general method and new technical tricks used in the 
derivation. Section II is devoted to the standard reformulation of the considered problem in terms of one-dimensional 
iV-particle system of quantum bosons with attractive ^-interactions @ . Here it is shown that the calculation of the 
the free energy probability distribution function, eq.Q, reduces to the summation over all the spectrum of eigenstates 
of this iV-particle problem. This summation is performed in Section III, where in the thermodynamic limit, t — > oo, 
the result, eas. (fTUl) - (|lU is derived. The concluding remarks and as well as the key points of the calculations are listed 
in the final Section IV. 

II. MAPPING TO QUANTUM BOSONS 

In terms of the partition function Z, eq.0, the function W(f), eq.©, can be defined as follows: 

W(f) = Km ~WT e M^Nf ) Z" (13) 

Ar=0 

where (...) denotes the averaging over quenched disorder. Indeed, substituting here eq.Q, we have 

00 I i\N r +°° 



w (f) = E -jJr / d f p (/') ex p{ AAr (/ ~ /')} (i4) 

/+oo 
df P(/')exp[-exp{A(/ -/')}] 

df P(f) 0(f - f) 



X- 



oo 
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which coincides with the definition, eq.©. 

Later on we will see that the integration over x in the definition of the partition function, eq.(J3]), requires proper 
regularization at both limits ±00. For that reason it is convenient to represent it in the form of two contributions: 

dx Z(x) + dx Z(x) = Z { _) + Z {+) (15) 

-00 Jo 

Thus, taking into account the definition eq.©, we get 

~ (-1)" r,*r, , B *T*.« tv . . , V A» 



W(f) = lim V l_J_ C xp{AJV/ + /3JV/of}(Z(-) + Z { 

A— »oo ' * 1\ ! 



(+)) 



N=0 
°° K.L=0 



00 (_l)K+L 

= } [m E J n exp{\(K + L)f + f3(K + L)f t}x 
\->oo * — ' K ! L\ 
K,L=0 

/0 />+oo 
dxi...dx K / dyi...dy L ^f(xi,...,x K ,y L ,...,y 1 ;t) 
-oo Jo 



(16) 



where 



*(xi,...,XAr;i) ee Z(.xi) Z(a; 2 ) ... Z(xat) (17) 
Using the relations, eqs.([T]), ([2]) and ((4]), after simple Gaussian averaging we obtain 

N T r 4> a (t) = X 

exp(-[3H N [(f)i,4>2, ...,4>n}) (18) 



a=l 



0a(O)=O 



*(aci, ...,x N ;t) = Yl 

where 

h n [<i>i,<i>2,-,<i>n] = o / ^(^[^^(t)] 2 -^^^^^-^^ 

^ - 70 \«=1 a^b 



rt / N N 



(19) 



The propagator <]/(x;i), eq. (fT8)) . describes iV trajectories ^(t) all starting at zero (<j> a (0) = 0), and coming to N 
different points {xi, ...,xn} at r = t. One can easily show that ^(x; t) can be obtained as the solution of the linear 
differential equation 



N N 

/?Wx;i) = - Y^dl^i^t) + - K J2t(xa-x b )n*\t) (20) 

a=l a^6 

with the initial condition 

*(x;0)=n^ =1 «( a ; o ) (21) 
and the interaction parameter k = /? 3 it. One can easily see that Eg. (1201) is the imaginary-time Schrodinger equation 

-^$(x;i) = ff«(xit) (22) 

with the Hamiltonian 

AT iV 

H = -^J2 d l-2 K ^ 5{Xa - Xb) (23) 

a— 1 a^b 

which describes TV bose-particles interacting via the attractive two-body potential —kS(x). A generic eigenstate of 
such system is characterized by N momenta {q a } (a = 1, ...,N) which are splitted into M (1 < M < N) "clusters" 
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described by continuous real momenta q a (a = 1, M) and having n a discrete imaginary "components" (for details 
see [HEMl): 



Qr = la - y("« +1 - 2r) 



(r = 1, ...,n a ) 



with the constraint 



A I 



(24) 



(25) 



Ot = l 



A generic solution ^(x, t) of the Schrodinger equation (|2"0)l with the initial conditions, Eq. (j2"Tj) , can be represented in 
the form of the linear combination of the eigenfunctions \l/q M ^(x): 



N 1 r r 



M=l 



where we have introduced the notation 



|C M (q,n)| 2 *W(x)*£O*(0) e xp{-£ M (q)*} 



r. M T /.-f co j 00 "I M 



(26) 



(27) 



and 8{k,m) is the Kronecker symbol; note that the presence of this Kronecker symbol in the above equation allows 
to extend the summations over n a 's to infinity. Here (non- normalized) eigenfunctions are [Tol |20| 



1 + IK 



sgn(o: a - gt) 

qv a - qv b 



exp 



(28) 



where the summation goes over N\ permutations V of N momenta q a , eg. (1241) . over iV particles £ a ; the normalization 
factor 



M I j K / \ 1 2 

ic^q, n)i» = — z— - n 



and the eigenvalues: 



E 



M 



l N x M 

(q) = ^E«« = »" 



24,3 ^ 

a— 1 a— 1 a — 1 



(29) 



(30) 



2/3 f-; a 2(3 

Note that the eigenfunctions, eq.(|28p. are symmetric with respect to permutations of all its arguments X\, ...,xn, and 

¥^\0)=M (31) 

In this way the problem of the calculation of the free energy probability distribution function, eq. (|16[) . reduces to 
the summation over all the spectrum of the eigenstates of the ./V-particle bosonic problem, which is parametrized by 
the set of both the continuous, {qi, qui}, and the discrete {n%, i%m}', {M = 1,...,N); (N — l,...,oo) degrees of 
freedom. 

III. FREE ENERGY PROBABILITY DISTRIBUTION FUNCTION 



Substituting eqs. (l2l)l) - (|3"Tj) into eq. (fT6"|) . (defining fo — jrf3 4 u 2 , the factor /o drops out of the further calculations) 



we get: 



oo 



W(f) = 1+ lim Yl 



K+L pA(K+L)/ 



K+L 

>< E^n 

M=l a=l 



K+L>1 
At 



E 



dq a 

2-KKTlr 



M 

6(j2 n * , N) \Cm{i, n )\ 2 Ik,l{<1, n) 052) 

a=l 



5 



where 



|C M (q,n)| 5 



M 

n 

a</3 



q a -qp - % -j{n a + np)\ 



(33) 



and 



j*.*(q,n) = e ee nn 

■p(K,L) -p(K) -p(L) a = l C=l 



A" 



K 

n 

a<6 



n 



dxi...dx K expUE(gpW - «e)^a 

oo<xi<...<2:k<0 a— 1 Q 



0<j/_L<...<yi<+oo 



dVL—dyi exp^E(^(i) + ie)j 



c=l 



(34) 



Here the summation over all permutations V of (K + L) momenta {qi, qic+L,} over K "negative" particles 
{xi, xk} and L "positive" particles {yL, ■ j/i} are divided into three parts: the permutations PW of K mo- 
menta (taken at random out of the total list {qi, qK+L,}) over K "negative" particles, the permutations of 
the remaining L momenta over L "positive" particles, and finally the permutations V^ K,L ' (or the exchange) of the 
momenta between the group "K" and the group "L". Note also that the integrations both over x a 's and over y c 's 
in ea. (|34p require proper regularization at — oo and +oo correspondingly. This is done in the standard way by intro- 
ducing a supplementary parameter e which will be set to zero in final results. The result of the integrations can be 
represented as follows: 



i K M*n) = e nn 

■p(K,L) a = l c=l 



E 



q v {K) - q p (L) - ik 



K 

n 



r„(-) 



q v (K)) a <b L 



,(-) 



q v (K) q v (K) 



p(') i/pCtH'tpCt) T^jijj-^ji) -r ... 

where 

?i ±} = <7q ± «e 
Using the " magic" Bethc ansatz combinatorial identity [l3| , 



n 



rJ+) 

q v (L) 



q% + lK 



(7 (+) -0 (+) 



E 



i 



N 

P qpMvi +q P 2)-(q P i + - + J|, L 



9pa 9pt lK 



1 



JV 



(where the summation goes over all permutations P of N momenta {gi, ...,g/v}) we get: 

i 'v t z- \ - TT TT gp< g) ~ 1yj L) ~ lK 



q a + qb + iK 
q a + qb 



(35) 



(36) 



(37) 



IK.L 



(q, n ) 



e nn 

■p(K.L) o=l C=l 



n K o (_) - LJ - w~ J 

lla=l ^(A') a<fcl_ "WA) T v k 



n 



r „H 



(-) 



,(-) 



llc=l %p(l) c<d 



a (+) 



(38) 



Further simplification comes from one important property of the Bethe ansatz wave function, ea. (|28|l . It has such 
structure that for ordered particles positions (e.g. x\ < X2 < ■■■ < xn) in the summation over permutations 
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the momenta q a belonging to the same cluster also remain ordered. In other words, if we consider the momenta, 
cq. (f2"4"|) . of a cluster a, {q", q% , q" a }, belonging correspondingly to the particles {x i± < x i2 < ... < %i„ a }, the 
permutation of any two momenta and q" of this ordered set gives zero contribution. Thus, in order to perform 
the summation over the permutations p^ ,L ' in eg. (1381) it is sufficient to split the momenta of each cluster into two 
parts: {qf, Qm \\lm 9n Q }i where m a = 0, 1, n a and where the momenta q", belong to the particles 

of the sector "if", while the momenta q^ a+1 ... 7 q" a belong to the particles of the sector "L". 
Let us introduce the numbering of the momenta of the sector "L" in the reversed order: 

C -> 1 1 

a *ck 



1m a + l 



where m Q + s a = n Q and (s.f. eq. (|24[) ) 



q*r = Va + y (n a + 1 - 2r) = <7 Q + y (m Q + s Q + 1 - 2r) 
By definition, the integer parameters {m a } and {s Q } fulfill the global constrains 



(39) 
(40) 



A I 



m a = K 



^ Sq = L 



(41) 



(42) 



In this way the summation over permutations 7 ,< - /f,L - ) in cq. (|3"5)l is changed by the summations over the integer 
parameters {m a } and {s a }: 



A I 



e (-) - n 

■p(K.L) 



a=l 



E * 

m a -\-s a >1 



M M 

(5(E m a , if) 5 (E Sa ' L 



ce=l 



(43) 



which allows to lift the summations over K, L, and {n Q } in eq. (|32j) . In terms of the parameters {m a } and {s a } the 
product factors in eq. (|38|) are expressed as follows: 



K M m a 

a—l a—1 r—1 



L M s a 

(+) T7 TT *«(+) 



II'/ 



a=l 



nn*- 

a— 1 r—1 



(44) 
(45) 



n „(-) +!(-) 



n n 

a=l Kr<r' 



n nn 



of-1 1 

1<q</9 r=l r'=l L ' + <£' 



(46) 



n 

c<d L 



„(+) , „(+) 

(+) +"(+) 



AI 8c, 

n n 

a— 1 Kr<r' 



>«(+) 



g »a(+) +g *a(+) 



n nn 

l<a</3 r=l r'=l 



g»? (+) + g *£ (+) -i« 



(47) 



if L 

nn 



q v (K) — q v w — i>K 
q v {K) — q v {L) 



n nn 

l<a<8 \r=l r'=l 



9r 



3 r> 



q? + q*;> 



nn 

r=l r' = l 



IK 



qr - q' r > 



M m a 



nnn 

a— 1 r—1 r' — l 



q*r' 



(48) 
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Substituting eqs. (l43|) - (|48|) into eq. (|38|) . and then substituting the resulting expression into eq. (|32t we obtain 



W(f) 



M=l a=llm a +s a >l J 



dq a 



G(q a ,m a ,s a ) 
2itK,(m a + s a ) 



|CW(q,m + s)| 2 Y[ S ap (q,m,s) 

l<a</3 



(49) 



where 



|C A f(q,m + s)| = ]]_ — -± 

a <{3 \Qa -90- —\ m a + S a + TUp + Sp)\ 



(50) 



g? (-) + g°> (-) +i« 



n 

r<r f 



g*° (+) + g*°, (+) -i« 



nn 



g" — g r , — IK 



(51) 



and 

^=nn 

r=l r'=l 



nn 



. , n *a(+) I ( + ) 

r=l r'=l L g r T g r i 



nn 



?r ~ q*r' ~ iK 



r=lr'=lL Qr q r' J r=lr'=l 



*/3 



Sa m l3 

nn 



g ,. ~ g£ ~ IK 



q*? ~ € 



The product factors in eq. ([5T1) can be easily expressed it terms of the Gamma functions: 

I 5 



n*" ( 5 = n[ &H>_ Y(™ Q +s Q +i)+iKr] = (*«) 



r=l 



s a —?n a iq a 

2 



1 s a +m a iq a ^ ^ 



(52) 



(53) 



(+) 



— (m a + s a + l) — z^r 



'(' 



m Q —s a i ^gQ 



p / 1 m Q +s Q 

1 2 2 



(54) 



n 



2-(m a -l)_ 



r m a - s a - 



r l- 



m a +s 



,(-) 



r i 



2iq a (- 



(55) 



n 



*a(+) 



*a (+) 



( + ) 



= 2 -(. -l) 



r( Ja -^ + ^)r(i-^ + *gl) 



m Q + 



2iq a W 



(56) 



nn 



— q T , — ik 



rfi + m 



T(l + m a )T(l + s c 



r=l r' = l L 

Substituting the above expressions into eq.(|51j) and using the standard relations for the Gamma functions, 

T(z)T(l-z) 

T(l + z) 

r (H 



sin(7rz) 
zT(z) 

v^r(i + 2z) 
2 2z r(i + z) 



(57) 

(58) 
(59) 

(60) 



for the factor Q, eq. (|5Tj) . we get 

G(q a ,m a ,s a ) 



'(s a + ^q a { - ] ) r(m a - f g Q (+)) T(l + m Q + s a ) 



2(«*-+»<.)r(m Q + s« + f g a <->) r(™ Q + s a - f r(i + m a )r(l + s a ) 

Similar calculations for the factor Q a p yield the following expression 



<?Qfl(q,ni, s) 



1+ m.^-.«-., _^ fe( - )+?/)( - ) 



1 _ ma+ mf)+Sa+Sff _ k ( ga (-) + g/? (-) 



1 + 



2 



(+) 



771 a +77l£+,S a — S£ 

2 



777. a + m p — S a. + S p 



1+ +i(q a -q ) 



- m a + m p + s q- — s p 



+ i(q a -qp) 



1p) 



ma+mp+Sa—Sfi 
2 



q&) 



i , m. a +m,fj+s a +sp 
1 ' 2 



£(««-9/8) 



1 ~r 2 J™ 



1 + 



m a -\-mp — s a -\-sp 



i + 



Redefining 



with 

the normalization factor (q, m + s)| 2 , eq. ([50|) . can be represented as follows: 

, , l2 ¥r |A(m a + s a ) - Xims + sA - ip a + ipp? 

|C A /(q,m + s)|- = || ' - 



a<f} \\(m a + s a ) + \(mp + 

M 

J^J [2X(m a + s a )] x det 



iPa + IP/3 



a=l 



\{m a + s a ) - ip a + \(mp + sp) + ipp 



a,0=l,...,M 



where we have used the Cauchy double alternant identity 



= (_l)^(M-D/2 det 



a a - b f . 



a,/3=l,...M 



with a a — p a — i\(m a + s a ) and b a — p a + i\(mp + sp). 

After rescaling, eq. (p5|) . for the exponential factor in eq. p§)) we find 

t n 2 1 

- ^g( TO « + S a)la + ^g( m « + S «) 3 + K m a + S a )f = - \(m a + S a )p a + -\ 3 (m a + S a ) 3 + \(m a + S a )f 
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The cubic exponential term can be linearized using the Airy function relation 



exp 



^X 3 {m a + s a ) 3 



+ 00 



dy a Al(y a ) exp A(m a + s a ) Da 



(68) 



Substituting eqs. (j6"5|) . (|()T)) and (p5|) into eq. (|4"9")l , and redefining y a — > y a + p 2 a — /, we get 



W(f) = lim { 1 + e 



M=l 



(-1) M 

M! 



M 

n 

ce=l 



+ dy a dp a , 2 r\ 

-rv-i 



x £ (-ir" +s °- 1 exp{A(m Q + SQ )y ct }^(^, 

m a -\-s a > 1 

x deti^[(ATO Q , As Q , p Q ); (Am^j, As/j, p /3 )] a (3=lj M x J| m, s 



l<ct</3 



where 



AT [(Am, As, p); (Am', As', p')] = 



1 



Am + As — ip + Am' + As' + ip' 



(69) 



(70) 



The crucial point of the further calculations is the procedure of taking the thermodynamic limit A — > 00. First of 
all one can easily note that according to eas. (|6"3")) and (|6"9")l . it is the parameters p a ~ Xq a which remain finite in the 
limit A — > 00. In other words, all the parameters q a which are not multiplied by A (e.g. in the expressions for Q a and 
Qafi-, eqs (|61[) and ([62j) ) have to be taken to zero in this limit. Simultaneously, the summations over {m a } and {s a } 
have to be performed. The general algorithm of such summation is in the following. Let us consider the example of 
the sum of a general type: 



M 



R(y,p) = lim Yl 



(-1)"° 1 exp{Xn a y a } 



p, -, An; n 



(71) 



where $ is a function which depend both of An a 's and n a 's (which are not multiplied by A). The summations in the 
above example can be represented in terms of the integrals in the complex plane: 



M 



R(y,p) = lim ]J 



1 



2i J c sin(7T2; Q ) 



exp{Az Q y Q } 



p, -, Az; z 



(72) 



where the integration goes over the contour C shown in Fig. 1(a). Shifting the contour to the position C shown in 
Fig. 1(b) (assuming that there is no contribution from 00), and redefining z — > z/X, in the limit A — ¥ 00 we get: 



R(y,p) = II 



1 

2Tri 



dz 



-exp{z a y a } 



C z a 



lim <fr 

A— >oo 



P z 

P ' A ' Z; A 



(73) 



where the parameters y a , p a and z a remain finite in the limit A — > 00. 

Let us consider now the summations in eq. (|69|) . Here the double sum can be represented as follows: 



E (-l) m - + "'- 1 /(m; B ) = E i- 1 )™*- 1 f(™;s)\ Sa=0 + E(-l)"- _1 /(m;8) 



m a =0 



m a -\-s a >1 



E (-lr^Ec-^^/^s) 



(74) 



m a —l s a — 1 

Thus, according to the above summation algorithm, we get 



lim E {-l) ma+s °>- 1 /(m;s) = 



1 



A— »oo 



m a -\-s a > 1 



(27Ti) 2 



(2m)S(z 2a ) + (2m)S(z la ) - l] Urn (75) 
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c 



1 2 3 4 5 



1 2 3 4 5 



(a) 



(b) 



FIG. 1: The contours of integration in the complex plane used for summing the series: (a) the original contour C; (b) the 
deformed contour C'\ 



where the rescaled integration parameters z\ a and z 2a remain finite in the limit A — > oo. Finally, taking into account 
the Gamma function properties, r(z)|| z |^ = l/z and T(l + z)l|zKo = 1) f° r the factors Q and Gap, eas.(f6"T j) - (j62|) , we 
easily find 

Um G(q a ,m a ,s a ) 2^+^ = lim q{% ^, ^W^-A+WA) 

A— >oo A— too \ A A A / 

(Zla + z 2a + iPa~ } ) (zi a + Z 2a ~ ipi +) ) 



and 



where 



lim 

A— ^oo 



[z 2a + ipi ') (z la - 

g Av m ' s ) = ItJAv T' t 



G* (Pa,Z la ,Z 2c 



P 



(76) 



(77) 



(78) 



Thus, in the limit A — >• oo the expression for the probability distribution function, eq. (|69[) . takes the form of the 
Fredholm determinant 



w(f) = i+E 



(-i) 



M 



..v n 

M=l a=l 



dy a dp a 
2ir 



Ai(y Q +Pa-/) 



(27T«) 



(27tz)«5(z 2q ) + (27rz)«5(z lQ ) - 1 



Z2a + iPa ' 



^2 a 



_ • (+) 



X det AT[(zi Q ,2;2a,Pa); (^1/3, ^2/3,P/3)] ct / 3 =1 



0=1,. ..,Af 



= det[l - K] 
with the kernel 



k[(zi, z 2 , p); (zx, z 2 ', p')] 



Z\ + z 2 — ip + z\' + z 2 ' + ip 



In the exponential representation of this determinant we get 



W(f) = exp[- ]T jj Tr /• 

M=l 



M 



(79) 



(80) 



(81) 
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where 



TyK 



M 



M 

n 

a=l 
1 



+ 00 



dy a dp a 
2ir 



(2-nif 



M 



dz la dz 2c 



& z l a z 2 a 



n 



a = l 



(2m)S(z 2a ) + (2m)5(z la ) - 1 



1 



zi< 



^2 Q + «pi } 



^2a 



_ ■ (+) 



z la + z 2a ~ *Pa + z la + l + z 2a+l + l Pa+l 

Here, by definition, it is assumed that Zi M+1 = Zi t (i = 1, 2) and pu+i = Pi- Substituting 
1 



(82) 



Zla + z 2a - iPa + z la+l + z 2a+l + iPa+1 

into ea.(|82p. we obtain 

M 



duj a exp 



- [Zl a + Z 2a - ip a + z l a +l + z 2a+l + iPa+l) Ua (83) 



TrK M = / duj 1 ...duj M TT 
J a=l 

where, by definition, luq = cum, and 
I f f dzi dz 2 



+00 



dy a dp a 
2tt 



Ai(Va+Pa + Ua + u a -i - f) exp{ip a (u! a - u a -i) } S(p a , y a ) 



S (p, y) = 



(2m) 



C z l z 2 



(2m)6(z 2 ) + (27r»)<y(*i) - 1 



Zl 



z 2 + ip 



z 2 



z\ — 



i(zi+z 2 )y Q 



Simple calculations yield: 



S(P,V) = ±[^[1 



Zl 



2ni J C i z\ V i(p — ie) 
I /' f dz\ dz 2 ( 



exp{ziy} + -!- / — (l - Z2 , ) cxp{z 2 y} 



Zl 



2iri J c , z 2 V i(p + ie) 



(2iti) 2 J J c , z\ z 2 \ z 2 + i(p — ie)/\ z\ — i(p + ie) 



1 



z 2 



exp{(zi + z 2 )y} 



(84) 



(85) 



1 



1 



-i(p — ie) i(p + ie) 
Taking the limit e — ¥ we find: 

Substituting this result into eq. (|84]) we obtain 



S(p, y) = S(y)5(p) 



poo M 

/ duji ... duj M TT Ai(w Q + w Q _i - /) 
Jo i L 



TtK m = / dui 



In other words, the free energy distribution function of our problem is given by the Fredholm determinant, 

W(f) = detfl - K- 

with the kernel 

K- f (w,w') = Ai(w + a/ - /) , (w, w' > 0) 

which is the GOE Tracy- Widom distribution [H, HI]. 



(86) 



(87) 



(88) 



(89) 



(90) 
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IV. CONCLUSIONS 

In this paper we have presented sufficiently simple derivation of the GOE Tracy- Widom distribution function for 
the free energy fluctuations in random directed polymers with free boundary conditions. The main message of this 
somewhat technical work is not the final result itself (which is not new anyway), but the demonstration of the efficiency 
of the general method and new technical tricks used in the derivation. By mapping the original problem to the N- 
particlc quantum boson system with attractive interactions the derivation is done in the framework of the integer 
replica series summations and the Bethe ansatz formalism for the quantum boson system. 

The key technical tricks of presented calculations includes the following points. First of all, to make the integration 
over particle coordinates of the Bethe ansatz propagator well defined one has to introduce proper regularization at 
±00 which requires formal splitting the partition function into two parts: the one in the positive particles coordinates 
sector (up to +00) and another one in the negative particles coordinates sector (down to —00), cqs.(fT5 |) - ([T6l) . Next is 
the "magic" Bethe ansatz combinatorial identity, eg. (1371) . which allows to perform the summation over the momenta 
permutations and "disentangle" sophisticated products containing in the Bethe ansatz propagator. One more trick is 
the reformulation of the summation over permutations of the momenta between the positive and the negative particles 
positions sectors in terms of the series summations, eq. (|43[) . which allows to represent the probability distribution 
function in terms of the problem of the series summations, eq . (|4"9")l . Finally, the crucial point of the considered 
derivation is the procedure of the series summations in the thermodynamic limit t — > 00. In this limit, due to the 
integral representation of the series, eas. (|7Tj) - flT3"l) . one obtains dramatic simplifications of some factors, eas. ([7p| - (l771) . 
in the expression for the probability distribution function, which allows to represent it in the form of the Fredholm 
determinant, eq. (j?9"1) . 

Hopefully, the experience obtained in the presented calculations would help to solve more serious long standing 
problems of this scope, such as the distribution function of the directed polymer's end point fluctuations or the 
statistical properties of the free energy fluctuations at different times. 
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